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Abstract. We show that any solution of a smooth first order dynamic equation can be made 
Lyapunov stable at will by the choice of an appropriate time-dependent Riemannian metric. 


1 Introduction 

The Lyapunov matrix of a first order dynamic equation is defined as the coefficient matrix 
of the variation equation, and its properties are used as a criteria of stability of solutions 
of this equation The problem is that the Lyapunov matrix is not a tensor under 

coordinate transformations, unless they are linear and time-independent. For instance, 
any Erst order dynamic equation dehnes a coordinate atlas such that its solution is con¬ 
stant on each coordinate chart (see Proposition |l| below) and, consequently, the Lyapunov 
matrix vanishes. 

We introduce the covariant Lyapunov tensor of a hrst order dynamic equation and 
study the Lyapunov stability of solutions when this tensor is negative-dehnite. The co¬ 
variant Lyapunov tensor essentially depends on the choice of a Riemannian metric. We 
show that, if a dynamic equation is defined by a complete vector held, the Lyapunov 
exponent of its solution can be made equal to any real number with respect to the ap¬ 
propriate time-dependent Riemannian metric. It follows that chaos in dynamical systems 
described by smooth (G°°) hrst order dynamic equations can be characterized in full by 
time-dependent Riemannian metrics. 

2 Geometry of first order dynamic equations in non-autonomous 
mechanics 

Let R be the time axis provided with the Cartesian coordinate t and transition functions 
t' = t-l-const. In geometric terms, a (smooth) hrst order dynamic equation in non- 
autonomous mechanics is dehned as a vector held 7 on a smooth hbre bundle 

( 1 ) 


1 


( 2 ) 


which obeys the condition yjdt = 1 |^, i.e., 

1 = dt + 7 ^( 9 a . 

The associated first order dynamic eqnation takes the form 

^ = 1 , ( 3 ) 

where are holonomic coordinates on TY. Its solntions are trajectories of the 

vector field 7 (|^). They assemble into a (regnlar) foliation T of Y. Eqnivalently, 7 (I) is 
defined as a connection on the fibre bnndle (|^. 

A fibre bnndle Y (|ID is trivial, bnt it admits different trivializations 

F = R X M, ( 4 ) 

distingnished by fibrations Y M. Moreover, Proposition below appeals to an atlas 
of F as a fibred manifold. If there is a trivialization (D snch that, with respect to the 
associated coordinates, the components 7^ of the connection 7 (j^) are independent of t, 
one says that 7 is a conservative hrst order dynamic eqnation on M. 

Proposition 1. Given a first order dynamic equation 7, there exists an atlas T = 
{{U]t, y“)} of a fibred manifold F —>■ R with time-independent transition functions y'°‘{y^) 
such that any solution s of j on each chart {U;t,y°‘) reads 

s°‘{t) = const., t G 7r{U) C R. 


Proof. The atlas T is an atlas of adapted coordinates for the foliation T of trajectories 
of the vector held (@) 0. □ 

Proposition 2. Let the vector field 7 be complete, i.e., there is a unique global 
solution of the dynamic equation 7 through each point ofY. Then there is a trivialization 
w of Y such that any solution s 0/ 7 reads 

s“(f) = const., t e R, 

with respect to associated bundle coordinates it,y°‘). 

Proof. If 7 is complete, the foliation tF of its trajectories is a hbration of F along 
these trajectories onto any hbre of F, e.g., F^^q — M. This hbration yields a desired 
trivialization |P, §] . □ 

One can think of the coordinates {t, in Proposition ^ as being the initial date 
coordinates becanse all points of the same trajectory diher from each other only in the 
temporal coordinate. 
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Let us consider the canonical lift 1^7 of the vector held 7 onto the vertical tangent 
bundle VY of 1 " —>■ M. With respect to the holonomic bundle coordinates on 

VY, it reads 

V"7 = 7 + dx = (5) 

This vector held obeys the condition V'j\dt = 1 , and dehnes the hrst order dynamic 
equation 


i = l, y^ ='y^{t,y''), 

Vt^ = di,-f^{t,y^)y^ 


(6a) 

(6b) 


on VY. The equation (pa]) coincides with the initial one (j|). The equation (pbD is the 
well-known variation equation. Substituting a solution s of the initial dynamic equation 
(|6a|) into (|6bf ), one obtains a linear dynamic equation whose solutions s are Jacobi helds of 
the solution s. In particular, if y —M is a vector bundle, there is the canonical splitting 
VY = Y xY and the morphism VY —> y so that s-|-s obeys the initial dynamic equation 
(fe) modulo the terms of order > 1 in s. 

Remark that, if 7 is the Hamilton equation for a Hamiltonian form H on 


y = V*Q ^ Q ^ M 


in time-dependent Hamiltonian mechanics, the variation equation 1/7 is also the Hamilton 
equation for the Hamiltonian form on V*VQ which is the canonical tangent lift of H onto 
VV*Q = V*VQ 


3 The covariant Lyapunov tensor 

The collection of coefficients 

= 5,7" ( 7 ) 

of the variation equation (§D is called the Laypunov matrix. Clearly, it is not a tensor 
under bundle coordinate transformations of the fibre bundle Y (p. Let us bring into 
a covariant tensor. 

The stability conditions of Lyapunov for solutions of a first order dynamic equation 

involve the notion of a distance between different solutions at an instant t. Therefore, 

let a fibre bundle y —M be provided with a Riemannian fibre metric g, defined as a 

2 

section of the symmetrized tensor product V V*Y Y oi the vertical cotangent bundle 
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V*Y of y —i> R. With respect to the holonomic coordinates (t, y^, y^) on V*Y, it takes 
the coordinate form 

9 = ^9ay{t, y^)dy^ V dy^, (8) 

where {dy^} are the holonomic hbre bases for V*Y. 

Given a hrst order differential equation 7, let 

G*7 = 7 - ( 9 ) 

be the canonical lift of the vector held 7 (|[) onto V*Y . It is a connection on V*Y —>■ R. 
Let us consider the Lie derivative of the Riemannian hbre metric g along the vector 
held 1/7 (H). It reads 

Lay = {Dtg)ay = Otg^y + j^dxgax + daj^gxy + Ogj^g^x- ( 10 ) 

This is a tensor with respect to any bundle coordinate transformation of the hbre bundle 
i). We agree to call it the covariant Lyapunov tensor. If g is an Euclidean metric, it 
comes to symmetrization 

Lay = < 9 „ 7 ^ + < 9 / 37 “ = 
of the Lyapunov matrix (|^). 

Let us point the following two properties of the covariant Lyapunov tensor. 

(i) Written with respect to the atlas T in Proposition the covariant Lyapunov tensor 

is 


Lab Ltgab' 


( 11 ) 


(ii) Given a solution s of the dynamic equation 7 and a solution s of the variation 
equation (|6b|) , we have 


Lan(t, s^(())s"s'’). 


( 12 ) 


The dehnition of the covariant Lyapunov tensor (pR]) depends on the choice of a Rie¬ 
mannian hbre metric on the hbre bundle Y. 


Proposition 3. If the vector field 7 is complete, there is a Riemannian fibre metric on 
Y such that the covariant Lyapunov tensor vanishes everywhere. 

Proof. Let us choose the atlas of the initial date coordinates in Proposition |^. Using the 
hbration Y —> Yt=o, one can provide Y with a time-independent Riemannian hbre metric 

9ab{t,y'") = h{t)g^^^{y'^) ( 13 ) 
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where Qabiy'^) is a Riemannian metric on the hbre and h{t) is a positive smooth 
function on M. The covariant Lyapnnov tensor with respect to the metric (]R) is 


Lab dth(]ab- 

Pntting h(t) = 1 , we obtain L = 0 . 


□ 


4 The local Lyapunov stability 

With the covariant Lyapnnov tensor (0)> we obtain the following variant of the stability 
condition of Lyapnnov. 

By the instantwise distance pt{s,s') between two solntions s and s' of the dynamic 
eqnation 7 at an instant t is meant the distance between the points s{t) and s'{t) in the 
Riemannian space {Yt,g{t)). Recall that a solntion s of a first order dynamic equation 
g is said to be locally Lyapunov stable at an instant f G M if, for any £ > 0 , there is 
h > 0 such that pt{s, s') < S implies pt'{s, s') < e for each t' which belongs to some half 
closed interval [f,.). Being locally Lyapunov stable with respect to g, a solution is so with 
respect to any Riemannian hbre metric on Y. 

Proposition 4. Let the covariant Lyapunov tensor L at a point y & Y be a negative- 
definite bilinear form on VyY. Then there is an open neighbourhood Uy of y and an open 
tubular neighbourhood Ug C Uy of the trajectory s through y in Uy such that 

Pt'>t{s,s') < pt{s,s') ( 14 ) 

for any t' G [f,) C TiiUs) and s' crossing Ug. 

Proof. Since the condition and the statement of Proposition ^ are coordinate-independent, 
let us choose a chart {U ; t, y°‘) of the bundle atlas 4 / in Proposition |l| which cover the point 
y. With respect to these coordinates, the covariant Lyapunov tensor takes the form (pd^. 
Restricted to U, the leaf s of the foliation JF is an imbedded submanifold, and has an 
open tubular neighbourhood W, i.e., any trajectory in U intersecting W does not leave 
W. There is an open neighbourhood Uy (Z U of the point y E Y where the Lyapunov 
tensor L (0) holds negative-dehnite. Put W = W H Uy. The foliation tF restricted to 
the tubular W' dehnes its hbration 

c : fP' ^ IP' n 


and the corresponding trivialization 

w'^T^iw') X (ip'np^(,)). 
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(15) 



There is an open neighbonrhood V C fl W of y in which can be provided with 
the normal coordinates (x“) defined by the Riemannian metric g in and centralized 
at y. Let us consider the open tubular Us = If is a subbundle of the trivial 

bundle (p! 5 |) endowed with the coordinates (t,x°‘), t E (,). Without a loss of generality, 
put t = 7 r{y) = 0 . With respect to these coordinates, the solution s in Ug reads s“(t) = 0 . 
Let 


= const., t E {, ) = 7 ^{Us), 


( 16 ) 


be another solution crossing Ug. The instantwise distance pt{s, s') between solutions s and 
s' is the distance between the points (t, 0 ) and {t,u) in the Riemannian space {Yt,g{t)). 
This distance does not exceed the length 


r 1 


Ptis,s') = 


1 1/2 


J gah{t,Tu'')u'^v!’dT 

. 0 


of the curve 

x°‘ = rn“, r E [ 0 , 1 ] 

in the Riemannian space {Yt,g{t)). At the same time, we have 
Pi=o(s,s') =-pt^Q{s,s') 

The temporal derivative of the function s') (13) reads 

1 

1 


3tPt(s,s') = 






( 17 ) 


( 18 ) 


( 19 ) 


Since the bilinear form dtgah = Lab is negative-definite at all points of the curve (p!8D, the 
derivative (|^ at all points f G (,) is also negative. Hence, we obtain 

Pt'>o(s, s') < Pt'>o(s, s') < Po(s, s') = po(s, s'). 


□ 

The inequality (|I^ shows that a solution of a hrst order differential equation which 
obeys the condition of Proposition ^ is locally Lyapunov stable at an instant t. Moreover, 
one can say that it is isometrically stable at f, and is so in some neighbourhood of t. Given 
another Riemannian fibre metric on Y, this solution remains locally Lyapunov stable, but 
need not be isometrically stable. 
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5 The asymptotic Lyapunov stability 


Let 7 (^) be a complete vector field. With a minor modification of the proof of Proposition 
one can state the following. 

Proposition 5. Let s be a solution of the first order dynamic equation 7. If there exists 
an open tubular neighbourhood Ug of the trajectory s where the covariant Lyapunov tensor 
o is negative-definite at all instants t exceeding some to, then 


lim - /9t(s, s')] < 0 

t'^00 


( 20 ) 


for any t > to and any solution s' (Z Ug. 

One can say that a solution s in Proposition ^ obeys the isometric asymptotic stability 
of Lyapunov. Of course, it is also asymptotically stable, but so is any solution of a smooth 
dynamic equation as follows. 

Proposition 6. Any solution of any first order dynamic equation defined by a complete 
vector field is Lyapunov stable at each instant of time and is asymptotically Lyapunov 
stable. 


Proof. These properties obviously hold with respect to the Riemannian fibre metric (|I 3 [) 
in Proposition ^ where h = 1 , and so do with respect to any Riemannian fibre metric on 
Y. □ 


It follows that a solution of a first order dynamic equation is Lyapunov unstable only 
if this equation is given by a non-complete vector field 7 (|^). In this case, there exists a 
local (but not necessarily global) Riemannian fibre metric on Y such that a solution s of 
7 is locally stable at a given instant t. 

One can improve Proposition ^ as follows. 

Proposition 7. Let X be a real number. Given a dynamic equation 7 defined by a 
complete vector field 7 ^ and its solution s, there is a Riemannian fibre metric on Y 
such that the Lyapunov spectrum reduces to A. 

Proof. Recall that the (upper) Lyapunov exponent of a solution s' with respect to a 
solution s is defined as the limit 

K{s, s') = lim ^ ln(pi(s, s')). ( 21 ) 

Let provide the fibre bundle Y with the Riemannian fibre metric (^) in Proposition § 
where h = exp(Af). A simple computation shows that the Laypunov exponent (| 2 T|) with 
respect to this metric is exactly A. □ 
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is negative, the solution s is said to be exponentially Lyapunov stable. If there exists at 
least one positive Lyapunov exponent, one speaks on chaos in a dynamical system P], 
Proposition ^ shows that chaos in smooth dynamical systems can be characterized in full 
by time-dependent Riemannian metrics. 
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